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TASK  1  -  SETA  DENSITY  ANALYSES 


!♦  Magnetic  Storm  Perturbations 

The  first  of  our  data  analysis  efforts  involved  the  SETA-1  data  covering  the  March  22 
-  April  10,  1979,  period.  Color  visualizations  of  the  latitude  vs.  time  development  of  the 
density  response  to  varying  geomagnetic  activity  were  developed.  A  number  of  color  plots 
have  been  produced  in  both  geomagnetic  and  geographic  coordinates,  for  both  density  and 
cross-track  wind.  These  figures  have  been  provided  to  Mr.  Marcos  at  Phillips  Laboratory. 
These  visualizations  have  enabled  us  to  evaluate  optimum  periods  to  run  the  CTIM  model 
so  that  the  SETA  densities  can  be  used  as  means  of  model  validation.  These  efforts  are 
described  further  in  the  MODEL  DEVELOPMENT  section. 

2.  High-Latitude  Structures;  SETA-1 

Recognizing  that  the  largest  data/model  discrepancies  occur  at  high  latitudes;  and 
also  that  the  existence  of  "density  cells"  (sometimes  consisting  of  depressions  in  density 
during  high  magnetic  activity  -  a  feature  not  predicted  by  any  empirical  models)  seems 
to  be  generally  accepted;  we  have  initiated  an  investigation  to  see  if  such  features  Ccin  be 
uncovered  in  the  SETA  data  in  a  statistical  sense.  We  have  begun  this  work  on  the  SETA- 
1  data  base,  since  both  densities  and  winds  are  available  during  this  period.  Figures  1-6 
summarize  these  results.  Note  that  density  and  wind  plots  are  available  for  both  day  (N. 
Hemis.  andS.  Hemis.)  and  night  (N.  Hemis.  only).  These  figures  are  shown  in  geomagnetic 
coordinates,  but  they  are  also  available  in  geographic  coordinates.  The  data  are  normalized 
to  45  degrees  latitude;  densities  are  normalized  to  unity  -  and  winds  represent  the  winds 
with  the  value  at  45  degrees  subtracted;  therefore,  the  densities  are  normalized  to  unity 
at  45  degrees,  and  the  winds  to  zero  at  these  latitudes.  The  purpose  of  this  normalization 
is  to  examine  only  the  latitude  structure  poleward  of  45°.  This  will  allow  future  studies  to 
mix  different  seasons  and  levels  of  solar  activity  with  minimum  contamination. 

These  figures  reveal  a  number  of  interesting  features  which  are  not  completely  under¬ 
stood.  In  Figure  1,  the  latitude  structure  is  distinctly  different  for  the  3  Kp  levels,  and  a 
'depression"  or  'bell"  seems  to  be  visible  for  intermediate  Kp  levels.  The  daytime  winds 
(Figure  2)  do  not  reveal  much  dependence  on  Kp.  In  Figure  3,  the  highest  Kp  curve  shows 
the  lowest  density  compared  to  45°,  again  suggestive  of  a  cell.  On  the  nightside,  the  winds 
(Figure  4)  increase  in  intensity  with  magnetic  activity;  daytime  densities  in  the  S.  Hemis. 
all  show  evidence  of  a  depression  at  high  latitudes  (Figure  5).  And,  in  contrast  to  the  N. 
Hemis.  wind  results,  the  daytime  winds  in  the  S.  Hemis.  (Figure  6)  show  a  distinct  inrease 
in  intensity  with  Kp.  It  is  our  hope  that  by  comparison  of  these  results  with  climatological 
results  from  CTIM,  we  will  begin  to  understand  the  physical  mechanisms  underlying  these 
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Density  Normalized  to  Den200  at  45  Deg 


Plot  of  Nighttime  Normalized  Density  in  Geomagnetic  Coords  —  Binned  by  Kp 


Figure  3 
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Plot  of  Nighttime  X-wind  in  Geomagnetic  Coords  —  Binned  by  Kp 


Density  Normalized  to  Den200  at  45  Deg 


Figure  5 
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Plot  of  Daytime  X-wind  in  Geomagnetic  Coords  —  Binned  by  Kp 


Figure  6 
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structures,  and  to  incorporate  them  into  future  predictive  capabilities.  This  has  indeed 
been  performed  with  the  more  extensive  SETA-2  and  SETA-3  data  sets,  as  described  in 
the  following. 

3.  High-Latitude  Structures;  SETA-2  and  SETA-3 

Based  on  our  reasonable  success  delineating  significant  high-latitude  structures  using 
the  SETA-1  data  (March/ April,  1979),  and  considering  the  potential  contribution  of  high- 
latitude  structures  to  the  15%  barrier  in  model/data  discrepancies,  this  area  has  been 
further  pursued  through  analyses  of  the  SETA-2  and  SETA-3  data  sets  during  1982  and 
1983,  respectively.  For  ALL  of  the  data  available  for  these  two  missions,  we  have  binned 
the  densities  into  seasonal  and  Kp  bins,  as  indicated  in  the  attached  table.  The  seasons  are 
summer,  winter,  equinox,  defined  respectively  according  to  the  months  JJA,  NDJ,  MASO 
for  the  N.  Hemisphere,  and  with  a  6-month  shift  for  the  S.  Hemisphere.  Note  that  there 
are  on  the  order  of  600-800  orbits  averaged  for  low  magnetic  activity,  300-500  for  moderate 
activity,  and  roughly  100-200  orbits  averaged  for  active  magnetic  conditions. 

Figures  7,  8,  9  delineate  these  results  for  the  SETA  data;  Figures  10  and  11  for 
the  corresponding  MSIS  model  results,  where  point-by-point  calculations  and  binning  are 
identical  to  those  for  the  SETA  data  points;  and  Figures  12  and  13  provide  climatological 
results  from  the  Jacchia  1970  (J70)  model  corresponding  to  July  and  December,  1983 
conditions.  Note  that  all  densities  are  normalized  to  a  value  of  1.0  at  45  degrees  latitude; 
this  was  done  so  that  density  structures  from  various  solar  activity  levels  could  be  averaged 
meaningfully.  We  note  that  the  SETA  data  reveal  some  clear  latitude  and  Kp  dependences 
that  the  MSIS  and  J70  models  do  not.  These  unmodeled  variations  could  be  potentially 
important  for  future  density  forecasting  efforts. 

We  have  begun  to  investigate  the  sources  of  the  above  variations  by  running  the  CTIM 
model  for  similar  conditions,  and  preparing  the  same  plots.  These  results  are  depicted  in 
Figures  14  and  15  for  June  and  December  CTIM  runs  for  similar  solar  conditions,  and  over 
the  same  range  of  geomagnetic  conditions.  We  note  that  in  many  cases  the  sense  of  the 
latitude  variation  and  its  dependence  on  Kp  are  similar  to  those  exhibited  by  the  SETA 
data,  although  less  pronounced.  The  presence  of  these  variations  allows  us  to  examine  the 
CTIM  model  to  uncover  and  physical  mechanisms  underlying  the  effects;  moreover,  once 
this  is  done,  it  may  then  be  possible  to  "correct"  CTIM  to  better  reproduce  the  density 
measurements  —  this  type  of  process  is  at  the  heart  of  our  validation  efforts. 

One  must  be  cognizant  that  some  part  of  the  SETA  variations  is  due  to  in-track  wind 
effects  on  the  measured  accelerations,  which  cannot  be  separated  from  the  density-induced 
accelerations.  The  CTIM  model  has  been  used  to  provide  a  check  on  this  issue.  Figures  16 
and  17  are  the  same  as  Figures  14  and  15,  except  that  the  in-track  projection  of  the  CTIM 
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Table  1 


Kp<=3 

3<  Kp  <5 

Kp>=5 

Season 

Winter 

705 

438 

110 

Transition 

814 

494 

213 

Summer 

573 

264 

138 

Winter:  Nov.,  1982,  Jan.,Feb.,  Dec.,  1983 

Transition:  Sep.,  Oct.,  1982,  Sep.,  Oct.,  Mar.,  1983 
Summer:  Jun.,JuI.,1982,  Jul.,  Aug.,  1983 
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Figure  8 
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Figure  9 
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Figure  1 1 
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Figure  13 
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Figure  15 
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Figure  16 
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CTIM,Dec  15(with  winds) 
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Figure. 17 
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wind  vectors  have  been  included  in  the  depicted  densities.  It  is  noted  that  there  is  some 
small  tendency  for  the  winds  to  accentuate  the  structures  noted  above,  but  that  these  are 
relatively  small.  This  suggests  that  the  densities  inferred  from  SETA,  and  their  latitudinal 
and  Kp  variations,  are  correct.  However,  we  plan  to  perform  a  second-level  check  on  CTIM 
winds  during  the  SETA-1  interval  to  ascertain  whether  we  might  have  underestimated  the 
wind  correction  effect. 

MODEL  DEVELOPMENT 

1.  Validation  Process  for  Atmospheric  Density  and  Ionospheric  Change 

During  the  first  year  we  have  focused  on  a  four  day  period  (March  20-23, 1979)  during 
which  two  geomagnetic  disturbances  occurred.  The  goal  was  to  determine  the  cause  of  the 
plateau  in  the  accuracy  of  model  prediction  of  atmospheric  density  for  the  last  25  years. 
The  variance  between  model  and  atmospheric  density  data  has  remained  the  same  over 
this  time  inspite  of  the  increase  in  understanding  and  modelling  of  upper  atmosphere  pro¬ 
cesses.  We  have  simulated  this  four-day  period  forcing  the  model  in  the  simplest  possible 
way:  using  statistical  patterns  of  auroral  precipitation  and  convective  electric  field.  The 
magnitude  of  the  statistical  patterns  were  defined  by  the  auroral  hemispheric  power  in¬ 
dex  derived  from  TIROS/NOAA  observation  of  auroral  particle  precipitation.  This  is  our 
standard  automated  process  for  performing  a  simulation  of  a  particular  interval,  providing 
the  power  index  is  available. 

The  statistical  analysis  revealed,  as  we  expected,  a  similar  variance  for  the  model/ data 
comparison  as  can  be  found  by  using  standard  empirical  models  such  as  MSIS  and  Jacchia. 
There  are  two  reasons  for  the  results.  First,  the  magnitude  of  the  driving  force  is  not 
correctly  specified,  and  second  the  particular  time  sequence  on  the  pulses  of  energy  input 
is  not  timed  correctly. 

The  first  problem  can  be  seen  by  examining  a  series  of  orbits  following  the  first  injec¬ 
tion  of  energy.  The  data  clearly  shows  the  orbit  passing  through  a  density  hole,  possibly 
driven  by  divergence  of  the  neutral  dynamics  rather  than  heat  input,  but  this  is  yet  to  be 
verified.  Using  statistical  magnitudes  for  the  magnetospheric  sources,  the  depth  of  these 
density  holes  are  too  shallow.  Experience  with  alternative  forcing  functions  have  shown 
that  a  much  stronger  electric  field  magnitude  is  required  to  produce  the  observed  depth  of 
the  holes. 

The  second  problem  is  related  to  the  first.  Examination  of  another  sequences  of  orbits 
from  the  SETA  data  reveal  clear  propagation  of  energy  equatorward  following  an  energy 
injection.  The  sequence  is  repeatable  in  the  data  with  first,  a  density  increase  seen  at 
high  latitudes,  followed  by  a  density  bulge  at  midlatitude  and  low  latitude  on  subsequent 
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orbits.  The  propagation  is  consistent  with  large  scale  internal  gravity  waves.  In  the 
numerical  simulation  the  same  wave  features  are  evident  but  with  reduced  magnitude  and 
with  incorrect  timing  sequence. 

The  increase  in  variance  in  the  data  during  these  two  types  of  density  perturbations 
is  clearly  seen  in  the  data.  It  is  clear  that  the  increase  in  variance  is  caused  partly  by 
these  two  distinct  processes,  but  it  is  not  clear  if  they  are  the  only  processes  that  cause 
the  problem.  To  address  this  question  the  next  step  is  to  try  to  match  the  magnitude  and 
time  sequence  of  the  magnetospheric  drivers  to  the  observations.  This  task  is  currently  in 
progress.  The  other  cause  for  the  model/data  variance  is  the  spatial  distribution  of  the 
driving  sources.  If  a  density  hole  can  be  simulated,  but  is  not  in  the  correct  location,  it 
indicates  a  correct  magnitude  of  the  driving  electric  field  and  momentum  source  but  an 
incorrect  spatial  distribution.  These  questions  will  be  addressed  in  the  following  year. 

An  alternative  approach  is  to  validate  the  model  by  comparing  with  the  ionospheric 
response.  lonosonde  data  for  this  period  was  analyzed  to  quantify  the  magnitude  of  the 
midlatitude  negative  phase  and  its  regional  dependence.  This  dataset  is  currently  under 
scrutiny  to  determine  if  it  can  be  used  to  further  quantify  the  model  input  fields. 

2.  Semi-Annual  Density  Variation 

On  the  longer  time-scales  one  of  the  largest  differences  between  empirical  models  and 
density  data  is  the  magnitude  of  the  semiannual  variation.  This  variation  is  imposed  within 
the  empirical  models  although  no  mechanism  had  been  accepted. 

One  of  the  focii  for  the  model  development  during  the  first  year  has  been  analysis  of  the 
semi-annual  variations  of  atmospheric  density  as  simulated  by  the  physically-based  model. 
Observations  have  demonstrated  a  significantly  larger  density  at  equinox  as  measured  at 
a  fixed  satellite  altitude.  Only  part  of  this  increase  in  amplitude  can  be  attributed  to 
thermal  expansion  from  higher  globally-averaged  temperature  at  equinox.  A  requirement 
of  the  model  is  to  be  able  to  reproduce  the  semi-annual  variation.  If  such  a  variation  is 
not  simulated  by  the  model  as  a  natural  phenomena,  then  it  will  be  necessary  to  impose 
such  a  variation  artifically.  No  physically-based  model  has  been  able  to  demonstate  such 
a  variation. 

Analysis  of  solstice  and  equinox  model  simulations  under  identical  solar  flux  reveals 
a  substantial  semiannual  variation,  with  an  amplitude  similar  to  observations.  Part  of  the 
change  results  from  greater  globally- averaged  heating  at  high  latitudes  from  Joule  heating. 
The  implication  is  that  the  sum  of  Joule  heating  in  the  summer  and  winter  hemisphere  is 
less  than  the  sum  of  the  north  and  south  hemispheres  at  equinox.  The  greater  globally- 
averaged  Joule  heating  at  equinox,  and  the  thermal  expansion,  accounts  for  approximately 
one-half  the  modelled  semi-annual  vaxiation. 
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The  remaining  semi-annual  density  variation  suggests  an  alternative  mechanism.  It 
is  suggested  that  the  global  scale,  interhemispheric  circulation  at  solstice  mixes  the  major 
thermospheric  species,  analagous  to  a  huge  turbulent  eddy.  The  effect  causes  less  diffu¬ 
sive  separation  of  species  at  solstice  which  tends  to  raise  molecular  nitrogen  and  oxygen 
densities  and  reduce  atomic  oxygen  density,  compared  with  equinox.  The  increased  mean 
mass,  at  solstice,  reduces  the  density  scale  height  and  leads  to  less  atmospheric  density  at 
a  given  altitude.  It  is  suggested  that  this  compression  of  the  atmosphere  at  solstice  can 
explain  a  large  fraction  of  the  semi-annual  density  anomaly. 

Fmrther  analysis  reveals  that  the  vertical  profile  of  the  amplitude  of  the  semi-annual 
variation  from  the  mixing  mechanism  in  the  physically-  based  model  agrees  well  with  MSIS. 
The  implication  is  that  little  of  the  mean  level  of  the  semi-annual  density  variation  is  caused 
by  temperature  changes.  Much  of  the  variability  from  year  to  year  in  the  amphtude  could 
arise  either  from  different  globally-averaged  heating  rates  in  equinox  and  solstice,  or  from 
changing  amphtudes  of  heating  or  turbulent  mixing  from  the  lower  atmosphere  sources. 
Further  analysis  on  the  solar  cycle  trend  in  the  semi-annual  variation  will  also  be  pursued. 

OTHER  ACTIVITIES 

A  paper  entitled  “Thermosphere  and  Ionosphere  Dynamics  during  High  Solar  Activity 
Equinox  Conditions*  by  F.A.  Marcos,  T.  Fuller-  Rowell,  J.M.  Forbes,  and  M.  Codrescu 
was  presented  at  the  COSPAR  Assembly  on  July  15,  1996,  by  T.  Fuller-Rowell.  This  work 
primarily  involved  a  comparison  between  SETA  data  and  CTIM  predictions  during  the 
March,  1979  storm  period,  and  related  interpretations.  A  number  of  figures  depicting  the 
latitude  vs.  time  response  of  the  thermospheric  density,  as  revealed  by  the  SETA  data  and 
as  modelled  by  CTIM,  were  constructed  for  this  presentation. 
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